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1. Introduction 



Let S 1 = [— 7r, 7r] with opposite ends identified, fix < rg < r% < tt and define an annulus 
L = {(r,8) : r < r <n,0 G S 1 }. 



The functions M± 



given by 



M±(r,6) = ±(rcos6> - l,rsin0) 

map L into the plane, the images A± = M±(L) being centred at (—1, 0) and at (1, 0) respectively. 
Let (u,v) denote the usual Cartesian coordinates in R 2 . We assume that ro,r% are such that 
the annuli intersect in two disjoint regions, and denote the intersection region in which the v 
coordinate is positive by S+ and the other by £_. See Figure [1] 



— r 




-1 1 

Figure 1. The manifold A CM. 2 (shaded). 



U 



We denote A = A + U A_ and £ = S + U £_ 
functions are given by 



Inverses M^ 1 : A 



x S 1 to the above 
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M± x = ( ±u) 2 +u 2 ,tan 



-l 



u± 1, 

Define a iwisf map A : i — * i by 

A(r, 0) = (r,0 + 2vr(r - r )/(n - r )), 

remarking that A leaves invariant the boundaries of L and otherwise rotates points about the origin 
by an angle that increases with the radial coordinate. The twist function r i— > 27r(r — Tq)/(ti — ro) 
has derivative c = 27r/(ri — r ) and is afHne. We define two twist maps Y : A — > A in the plane, 
given by 

*(«,„) = / ^ oAo ^M X(u,v)eA + 

\ {u,v) otherwise, 



(M^oA-^Mz\u,v) if («,«)€ A- 
[ (u, v) otherwise. 

Definition (Linked-twist map). ^4 linked-twist map : A — > A is i/ie composition O — Y o 

We call this a linked-twist map m i/ie plane or say that it is planar. It preserves the Lebesgue 
measure \i on A. The relative direction of the two rotations (here they are opposite) affects the 
ergodic properties of linked-twist maps. Sturman et al. (2006) discuss this in some detail; in their 
terminology the present map is co-twisting. We illustrate in Figured 




Figure 2. One iteration of the planar linked-twist map. Part (a) shows some initial conditions 
in the form of a red horizontal line across the left-hand annulus A+ . Part (b) shows the image 
of these points under the twist map $ and part (c) shows the image under the linked-twist map 
6 = To <S>. 



The purpose of this paper is to prove the following: 

Theorem 1.1. Let r$ — 2 and r% = \fl . The planar linked-twist map : A — > A is isomorphic 
to a Bernoulli shift. 

We make some remarks. Crucial progress toward this result was made by Wojtkowski (1980) 
and our work builds upon his. He proved that the system considered here is amongst a family of 
such systems that possess an ergodic partition, using the technique of finding an invariant cone- 
field; we review this result in Section [3] The work lead him to conjecture that such systems also 
are mixing and our result shows that this is indeed the case for the system considered. We discuss 
the reasons for restricting to this one example in Section [8] 

Our paper is organised as follows. We discuss the recent resurgence of interest in linked- 
twist maps in Section [5] and the work of Wojtkowski in Section [3J The cornerstone of our proof of 
Theorem II .11 is the introduction of new coordinates for the manifold A and we do this in Section |4j 
Correspondingly we give a new expression for the planar linked-twist map in Section [5] In 
Section [S] we introduce a new invariant cone-field and show that it is preserved by the differential 
DQ. Unlike the cone-field introduced by Wojtkowski, ours affords us sufficient control over the 
orientation of local invariant manifolds to deduce strong ergodic properties; we give the details 
in Section [3 appealing to the work of Katok et al. (1986) to complete the proof. We make some 
concluding remarks in Section [8] 
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2. Background to the problem 

The study of planar linked-twist maps was motivated by a number of authors. Bowen (1978) 
showed that certain such maps have positive topological entropy, and asked whether they possessed 
any ergodic properties. Similar maps were shown by Braun (1981) to arise as an approximate model 
of the global flow for the Stormer problem, and were encountered by Thurston (1988) in his study 
of diffcomorphisms of surfaces. 

Considering briefly a more general linked-twist map Q j k = T k o for inte gers j, k (where 
j = k = 1 corresponds to the present case), Devaney (1978) showed that if jk ^ then there is an 
invariant, zero-measure Cantor set on which Qj t k is topologically conjugate to a subshift of finite 
type. Wojtkowski (1980) showed that under the same hypothesis, there are restrictions on the size 
of the annuli which guarantee that 0^ has an ergodic partition. The restrictions are stronger for 
the case jk < than for the case jk > 0; we give details for the latter case in Section [3l 

In an unpublished note Przytycki (1981) considers a variety of linked-twist maps including 
the present kind. In Przytycki (1986) he shows that under certain conditions periodic saddles and 
homoclinic points are dense for this large class of maps and moreover that they are topologically 
transitive. 

In recent years the study of linked-twist maps has taken on a new significance owing to 
developments in our understanding of the mechanisms underlying good mixing of fluids. Ottino 
(1989) has shown that the single most important feature to incorporate in the design of any fluid 
mixing device is a 'crossing of streamlines', by which we mean that flow occurs periodically in two 
transversal directions. That linked-twist maps provide a suitable paradigm for this design process 
was highlighted in Ottino & Wiggins (2004) and has been discussed at much greater length in 
Wiggins & Ottino (2004) and Sturman et al. (2006). 

This renewed emphasis on linked-twist maps in applications serves to motivate the nature of 
our research on this subject. While ergodic theorists have developed a very powerful and general 
framework for understanding the nature of ergodic behaviour in general dynamical systems (e.g. 
see the work of Liverani & Wojtkowski (1995)) there are very few situations relevant to applications 
where it is shown that the hypotheses necessary to conclude the existence of a particular ergodic 
property are satisfied for that particular example. This is essential for applications, and it is 
precisely in the spirit of our results. 

The situation is very reminiscent of the development of applied dynamical systems theory 
in the 1970s. Whilst it was know that generically stable and unstable manifolds of hyperbolic 
periodic orbits intersected transversely, and that the transverse intersections give rise to nearby 
Smale horseshoes, showing that this situation occurred in examples of interest to applications 
required significant further work (and research along these lines for concrete applications continues 
to this day). An excellent example illustrating this point is the work of Devaney & Nitecki (1979) 
on showing the conditions under which the Henon map possessed an invariant set on which the 
dynamics was conjugate to a shift map (i.e. the map possessed a 'horseshoe'). In that work 
estimates specific to the Henon map had to be carried out to show that the map satisfied the 
Conley-Moser criteria for the existence of such an invariant set, as given in Moser (1973). 

3. Wojtkowski's results 

Here we describe Wojtkowski's (1980) criteria for the planar linked-twist map to have an ergodic 
partition, defined as follows: 

Definition (Ergodic partition). O is said to have an ergodic partition if and only if A can be 
partitioned into at most countably many positive measure, Q-invariant, pairwise- disjoint sets Ai 
on which the restriction of O is ergodic. Moreover we require that each ergodic component will be 
the union of finitely many Bernoulli components which are permuted by the map, i.e. each set Ai 
has the form Ai — [J™^ Ai.j where for each j the restriction of O n W to Ai_j is Bernoulli. 

Let w = (u, v) € £ and denote by a(w) £ (0, tt) the angle at which the segment connecting w 
to (—1,0) meets the segment connecting w to (1,0). Let 



A Bernoulli linked-twist map in the plane 



4 



cot aH 

V = sup , (3.1) 

wE s r(w) 

where r(w) denotes the Euclidean distance from w to (— 1, 0). Recall that c = 2tt/ (ri — ro) denotes 
the derivative of the twist functions. Wojtkowski proved the following: 

Theorem 3.1 (Wojtkowski (1980)). If 

c> 2r, (3.2) 

then the linked-twist map : A — > A has an ergodic partition. 

In the same paper he conjectures that under the assumptions of Theorem 13.11 then also 
has the -FT-property. This would, by the work of Chernov & Haskell (1996), imply that it has the 
Bernoulli property. 

We discuss the proof of Theorem 13.11 briefly. It is easily argued that /i-a.e. w £ A lands 
in E under iteration of O and moreover returns to E infinitely many times, for those points not 
satisfying this condition must be rigid rotations around one of the annuli, and must have rational 
angle of rotation, else their orbit would be dense and hit E. From the strict monotonicity of the 
twist function one infers that such points are contained within a set of measure zero. 

Consequently for a full-measure set of points we may talk of the return map to E, or just the 
return map as we shall usually abbreviate it. Following Wojtkowski we define the return map on 
Af+ J (E) C L (rather than on E itself) as the map E : Af+ X (E) -> M+ X (E) so that 

S = M~ 1 o 0* o A/+, 
where i is the smallest (strictly) positive integer for which 4 (M+(r, 8)) £ E. 

dr 

TuiL 

de- ^ 




Figure 3. The invariant expansive cone U C T^L is shown in the left-hand figure. In the 
right-hand figure is the image of the cone under the differential map D9j; (dark-shaded) with 
the original cone (light-shaded) included for comparison. Observe how the cone is mapped into 
itself and vectors within it are expanded. 

For ui = (r, 9) e M+ X (S) let ft = dr, ft = d0 give coordinates in the tangent space T U L and 
define the cone 

Wojtkowski establishes that U is invariant under, and expanded by, the derivative -D0£. We 
illustrate the situation in Figure [31 More precisely, define the cone field 

U+= |J U(r,9) 



and let ||-|| be the norm in Tr Tt g\L induced by the Riemannian metric, i.e. ||(ft,ft)|| = \J fi\ + r 2 /?^- 
We have the following: 

Proposition 3.2 (Wojtkowski (1980)). D®j;(U+) C E7+. Furthermore there is a constant A > 1, 
independent of (r,9) or [3, and for vectors (3 £ U+ we have \\DQ^j3\\ ^ A|j/3|j. 

A detailed proof may also be found in Sturman et al. (2006). To arrive at the ergodic 
partition one determines that /x-a.e. point returns to E not just infinitely many times but with 
positive frequency, combines this with Proposition 13.21 to deduce non-zero Lyapunov exponents 
for such points and appeals to the theorem of Katok et al. (1986), which extends results of Pesin 
(1977) to certain non-differentiable systems. 
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4. Definition of the new coordinates 



Recall that S 1 = [— it, n] with opposite ends identified and that we take ro = 2 and T\ = \fl . Let 
X = [2, Vf\ and -X = [-V7, -2]. 

We introduce the new coordinates in two stages, starting with the annulus A + . With reference 
to the left-hand part of Figure [4] A + is divided naturally into three components: that part which 
intersects the annulus A_ (i.e. the region S, which we have light-shaded) and the remaining 
connected components A\ (dark-shaded) and A° + (unshaded), which respectively lie 'inside' and 
'outside' of the annulus A_ (not shown in the figure). 

We will provide first the definition and second some discussion. Let i^:Ix [0, 7r] — > [0, 7r] be 
given by 



20 



cos 



-1 1 



V7- 



^Jr 2 - 4r cos 9 + 4 
(tt - Vf) (e - cos" 1 : 

(TT-oos-i^a; 



if M + (r,6) e 
if M + (r,6) e £ 

if M+(r,0) € A° 



and extend ^ to a function ip : X x S 1 — * S 1 by insisting that it be an odd function of 9, i.e. that 
ip(r, -9) = -ip(r, 9). Finally let * : X x S 1 -> X x S 1 be given by 

#(r,0) = (r,tp(r,6)). 

Definition (New coordinates (x,y) on A+). Given (u,v) G define (x(u,v),y(u,v)) = 

i£> oM+ x (u,v). 




The definitions of "0 and ^ are somewhat opaque so let us now motivate them. Consider a 
point (u,v) G S + and denote (r, 9) = Af^ 1 (u,w). Then 
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(x, y) = *(r, 6) = (r, \J r 2 - 4r cos 6 + 4) = fr, ^/(2 



— r cos I 



\2 i r 2 s j n 2 , 



and substituting (u,v) — M + (r,9) = (rcosd — l,rsin0) gives 



(x, y) = (^/{l + uY + v\ V(l-«) 2 + « 2 ) • 



So x(it, v) is the Euclidean distance from (it, v) to (—1, 0) whereas y(u, v) is the Euclidean distance 
from (it, v) to (1, 0); coordinates defined in this way are often called two-centre bipolar coordinates. 
If we instead take (u,v) £ £_ then y becomes negative but still \y\ gives the distance to (1,0). 

The remainder of the definition of ip (i.e. for (r, 9) £ M^ 1 (A + \T 1 )) ensures that $ : I x S 1 ^ 
IxS 1 homeomorphically. This condition alone does not uniquely extend ip and so the definition 
given is just one of many possibilities. 

We extend the new coordinates to all of A, making use of the rotation i : T 2 — > T 2 given by 

t{x,y) = (~y,x). 

Definition (New coordinates (x,y) on A_). Given (u,v) £ y4_\S_ define (x(u,v),y(u,v)) = 
tof o MZ (u, v). Given (u, v) £ S_ define (x(u, v), y(u, v)) = —u o o Afr 1 (u, w). 

Again some discussion is required, in particular as to why we have introduced a minus 
sign for the case (u,v) £ £_. First observe that MZ 1 : A- -> X x S 1 and i : T 2 -> T 2 
are both homeomorphisms. So without the minus sign it follows from the definition of ifj that 
l o o MZ 1 : A_ — » S 1 x X homeomorphically. Moreover for (it, u) €E S± we have 



tof o M_ 



; x (n,u) = (±V(i + «) 2 + « 2 , V(i-") 2 + « 2 ) 



This agrees with our previous definition of (x, y) on S + whereas on £_ the minus sign must be 
introduced. Figure [5] illustrates A in the original Cartesians (it, v) and in the new coordinates 

(i^jeS'xS 1 . 




Figure 5. The manifold A, illustrated in its native Cartesian coordinates and in the new 
coordinates (x,y) £ S 1 X S 1 . Notice that there are two distinct representations of S— shown in 
the right-hand figure. The bottom-right representation is the correct one, the other (top-left) 
one is shown for sake of completion. 



5. The map expressed in the new coordinates 

We begin by showing that the linked-twist map considered (i.e. with r = 2 and r\ = y/7) has an 
ergodic partition; in fact we show that a whole class of linked-twist maps, including this one, have 
that property. The proof is simplified by the new coordinates introduced in the previous section. 
Following this we express in the new coordinates. 

Lemma 5.1. Let 2 ^ ro < r± ^ \fl . Then condition K3. 2)) is satisfied. 
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Proof. By symmetry it is enough to show that the condition holds on S + . Moreover the condition 
is implied by 

7T 

sup cot a(w) < —= . (5.1) 

«>GE + ' v 7 — 2 

Let (x, y) = 9 o M7 (to) give w G £+ in the new coordinates. The angle a appears in a triangle 
in which its adjacent sides have lengths x and y, and the opposite side has length 2 (the Euclidean 
distance between the centres of the annuli). The law of cosines says that 

T; 2 _|_ y 2 — 4 

cosa= - . (5.2) 

2xy V ' 

The partial derivative of (|5.2p with respect to x is given by 

d (cQga)= 1 Z/ 2 ~4 
(9a; 2y 2x 2 y 

and, using x.y £ 1, we calculate that l/2y G [l/2\/7, 1/4] and (y 2 - 4)/2x 2 y G [0,3/16]. It is 
easily checked that \j2\fl > 3/16 and so the derivative is always positive. Consequently cosct 
is an increasing function of x and so a is a decreasing function of x. By symmetry a is also a 
decreasing function of y. Combining these facts with (|5.2|) we find that a G [cos _1 (5/7), 7r/3]. 
Recall that cot is positive and decreasing on (0,7r/2) so that 

_■■ 5 5\/6 7r 

sup cot a = cot ml a = cot cos — = < — = 

„,es+ 7 12 \/7-2 

and the proof is complete. □ 
We now express O : A — > A in the new coordinates. Let 

i?= {V o M-\A + )} U (tofoM: 1 ^-^-)} 

as shown in Figure [6][a). There is a one-to-one correspondence between points in R and points in 
A. Let F : R — > R denote the map $ : A — > A in the new coordinates. Recalling our definition of 
$ in (jl.lj) we have 



F(x,y) = * o M+ 1 o $ o M + o * 1 (.T,y) = | 



f oAor'fi,!/) if elxS 1 
(a;,y) otherwise. 

F is a homeomorphism of R. Also let 

R' = {*oM+ 1 (A+\E_)} U {lofoMI 1 ^)} 

as illustrated in Figure [SJb). Again, there is a one-to-one correspondence between points in R' 
and points in A. Let G : R' — ► i?' denote the map T : ^4 — * A in the new coordinates. Recalling 
our definition of T in (|1.2p we have 



G(x,y) = io$o A 1 o* 1 ot 1 (x,y) = 



loF-'or'fi.y) if (ac, y) G S 1 x I 
(a;, y) otherwise. 

G is a homeomorphism of i?'. 

To compose -F and G we require the natural bijections R <— > R' , equal to —id on T x —I and 
id otherwise, and its inverse R' i— > i?. The representation of O in the new coordinates is thus given 
by 

H = (R' h-> R) o G o (i? R') o F = (R' h-> R) o t o F" 1 o tT 1 o (i?' ^ R) o F. 
To simplify the above expression let O, f2 _1 : S 1 x S 1 — * S 1 x S 1 be given by 

±1 , i ±l if (x, y) G I x =pl 



n ±A (a:,») = 



t Tl otherwise 
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then 



H = n~ 1 o F- 1 oflo F. 



Let S C R be the image of the 'intersection region' E, i.e. S = (I x I) U (I x —X). For z E S 
we define the return map Hs : S —> S, analogous to the return map 0£ . 



1 - 



-7T - 



(a) 



□ 



R 



— 7T — 7"! — J"o 







7"! 7T 



7T 



-ro 

-7T 



(b) 



+ 



— w—ri — rg 







Figure 6. The manifolds R, R' C T 2 = S 1 xS 1 , in parts (a) and (b) respectively . Each is in one- 
to-one correspondence with A but S_ is represented differently in each. F is a homeomorphism 
of R and G is a homeomorphism of R' . 



6. A new invariant tangent cone 

In this section we study the derivative DH Z for z € R. We approach this by studying the derivative 
DF^ 1 for z 6 R, observing that this is the identity when z^IxS 1 , so the only interesting case 
is for zelxS 1 . In that case 

F ±x (x : y) ^^oA^oT^y) = (x, (x, y) ± c(x - r ))) 

where : X x S 1 — > S 1 is defined by (ac, ^(x, y)) = y. To simplify the expression write 

y±=y±(x,y) =ip~ 1 (x,y)±c(x-r ) and f±(x,y) = ip (x,y±) . 
If Di, D 2 denote the usual differential operators then the Jacobians of F ±x are given by 

np±l = ( 1 \ 

^ D 1 f ± {x,y) D 2 f±(x,y) J ' 
The derivatives of /± are given by 

£>i/±0, y) = y±) + A^O*, ir±) [^i^te, y) ± c] , (6.1) 
D 2 f±(x,y) = D 2 il){x,y±)D 2 il)- l {x,y). (6.2) 

Let 61 = dx, 62 = dy give coordinates in the tangent space T z (S 1 x S 1 ) = T^T 2 to a point 
z = (x, y) € R, and define the cones 

C{g) = {(b 1 ,b 2 )ib 1 b a >0} and C{z) = {(61, b 2 ) : 6162 < 0} . 
The cone C is illustrated in Figure [71 Define the cone fields 

c + = U and c - = U 

The remainder of this section is devoted to proving that DH preserves the cone field C + . 
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Figure 7. The invariant cone C C T Z R is shown in the left-hand figure. In the right-hand 
figure is the image of the cone under the differential map DH$. The fact that C is invariant 
under this differential is immediately implied by Proposition 16, ll In contrast to the situation 
in Proposition 13, 21 we do not claim that this cone is expanded by DHg, although it will follow 
from the results of Section [7] that this is true on average. 



Proposition 6.1. Let r = 2 and r\ — V7. If ze Rand we C(z) c T Z T 2 then 

DH z w G C{H{z)) C T H(z) T 2 . 

Proof. Some observations will simplify the task. First, from the definition of H it is enough to 
show the result holds for each of DF Z and D (Jl^ 1 o F^ 1 o Q) . Second, by the chain rule and the 
easy observation that Dttf 1 each map C+ into C_ and vice versa, it is enough to show that DF^ 1 
preserves C±. Third, as these derivatives are otherwise the identity, it is enough to consider only 
zelxS 1 . Finally fourth, let b 2 ) G T Z T 2 and define 



DF 



±i 



l>2 



h 

hD 1 f±(x,y)+b 2 D 2 f±(x ) y) 



We have 



^- = D 1 f ± (x,y) + ^D 2 f ± (x,y), 



so it is enough to show that for each z — (x, y) E X x S 1 we have 
±D 1 f±(x,y)>0 and D 2 f±(x,y) ^ 0. 



(6.3) 



We claim that 
Ditp G 



D 2 i/j G 



D^ 1 G 



4- 



and D 2 %l) 1 G 



It is elementary, using (|6.ip and (|6.2p . to check that in this case (|6.3j) is satisfied. Unfortunately 
proving the claim will require some extensive calculations. To keep the length of the proof within 
reasonable limits we prove only the first assertion; the others are proved similarly. For full details 
see Springham (2008). We further restrict to the case 9^0; the case 9^0 follows by symmetry. 

Now, Diijj(r, 9) takes three different forms corresponding to M + (r, 9) in each of A 1 ,, S + and 
A ,. We deal with each case in turn. 



(i) Let M + (r, 9) G A\ . We will estimate the range of 

d_( 29 

Or 



cos 



2(c 



(6.4) 



1 16 



The numerator attains a minimum of zero and a maximum of 



sup I 

M + (r,6)eA l + 



sup COS 

rel 



cos 1 inf - 

r€X 4 



7T 

3' 



Observe that 
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• r -1 r -1 r -1 V7 5 

mi cos - = cos sup - = cos — — > - 



so the range of (cos 1 |) 2 falls within f§§,Trj- Notice also that ^ € [|> §) an< ^ so 
Combining these gives a range for (|6.4p of 0, p^Tf) c [P> f) - 



(ii) Let M + (r, 9) E S. We will estimate the range of 



\/2 _\/3 
2 ' 2 



9r 



4r cos 6 1 + 4 



2cos6» 



VT 



4r cos + 4 



(6.5) 



By design the denominator has range T. For the numerator observe that the angle 9 occurs 
in a triangle where the adjacent sides have lengths 2 and r and where the opposite side, call 
it r'j has length in I. Using the law of cosines we have 



r — 2 cos 9 = r — 
The partial derivatives 



■ 4 - r'' 



r 2 -4- 



2r 



2r 



8 



(r — 2 cos ( 



r'2 , d , 

and — — ( r — 2 cos ( 



r 
r 



dr y 2r 2 dr 

are both positive for r,r' £l and so the numerator is an increasing function of each. This 



gives a range for the numerator of 



1 5^7 

' 7 



and for the quotient (|6.5p of 



7 ' 14 



c(o,i). 



(iii) Finally let M + (r, 6>) g . We estimate the range of 




(r 2 +3) (7T-V7)(7T 



r(rr-cos-i^) ^- (r 2 - 9) (r 2 - 1) 



(6.6) 



The numerator is non-negative but may be zero when 9 — %, this giving a lower bound. An 
upper bound requires the calculation 



inf 1 

M + (r,6»)eA° 



= inf cos 

rex 



_ x r 2 - 3 
4r 



_i r 2 - 3 
cos sup 

rex 4r 



cos 



and one can check that n — 9 < 2 as a consequence. Observing that tt — V7 < | gives 
upper bound for the numerator of 10. For the denominator we need also the calculation 



an 



Then tt 



sup cos 

rel 



COS^ 1 G 



4r 



3 _ x . r 2 

= cos mf 



rex 4r 



4r 



7T — COS 1 i,7T — COS 1 



c 



1 



3.2) 



Simple calculus gives 



y/— (r 2 — 9) (r 2 — 1) e [2\/3, 4] for r G X and so the denominator takes values in (6\/3, 8-\/7)- 
Consequently the derivative (|6.6p takes values in [0, 5V3/9) C [0, 1). 

Collectively the three calculations show that D\i\) E [0, |), proving the first of the four claims. A 
similar approach yields the remaining parts and concludes the proof. □ 



7. The Bernoulli property 

We now conclude the proof of our main result. For ^-a.e. w E A Sturman et al. (2006) show how 
the work of Wojtkowski (1980) gives a positive Lyapunov exponent associated to w. The theorem 
of Katok et al. (1986) implies that there is an unstable manifold ^(w) c A and an unstable 
subspace E u (w) c T lu R 2 . Analogously there is a stable manifold and subspace. Moreover Katok 
et al. (1986) give the following condition as sufficient for the Bernoulli property: for a.e. w, w' E A 
and for all sufficiently large natural numbers m and n 
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e m (7"H)ne-"(7 s (i»'))^. (7.1) 

Let z, z' denote w, w' in the new coordinates, so z is either o M+ l (w) or t o f o MZ 1 (w) 
as appropriate and similarly for z' . For sake of discussion take z — \& o AI^ 1 (w) and define 7" (z) 
to be the maximal connected, smooth component of $oM7 1 (7"(m)) containing z. The piecewise 
smoothness of \1/ ensures that there are only finitely many such smooth components and so the 
one containing z will have positive length. By 'a.e. z £ R' we mean those z corresponding to some 
full /^-measure set in A. We prove the following which implies (|7.ip : for a.e. z,z' £ R and for all 
sufficiently large natural numbers m and n 

r'(7 a ( Z ))n_ff-"(7 s (z'))^. (7.2) 

There are two facts upon which our proof relics. The first is that the length, naturally defined, 
of H m (7" (z)) grows arbitrarily large with m. (Essentially this follows from Proposition 13.21 and 
from the one-dimensional mean- value theorem, although some care needs to be taken as H is only 
piecewise smooth. The full proof is omitted for reasons of length, but can be found in Springham 
(2008).) The second concerns the orientation of j u (z): 

Proposition 7.1. For a.e. z £ R we have E u (z) C C C T Z T 2 . 

Proof. Fix z at which there is a positive and a negative Lyapunov exponent and consider the 
tangent space Tjjn/^H 2 for some n £ N. This may be written as a direct sum E™ E^ where 
E% = E u (H n (z)) and E^ = E s (H n (z)) are one-dimensional unstable and stable subspaces 
respectively. Fix some vq £ C\Eq then by the invariance of stable directions and of C 
(Proposition 16. ip one has v n = DH n v £ C\J5*. 

We show that if n is sufficiently large then E™ C C which gives the result. Assume for a 
contradiction that one may find arbitrarily large n £ N so that the inclusion does not hold. Then 
invariance of C and of the unstable subspace imply that the it does not hold for any n £ N. Let 

v n = a n e% + b n e n 

where en U ^ £ En is an (un)stable unit vector and where a n ,b n are non-zero and, without loss 



of generality, positive. Uniform expansion of unstable vectors by the return map (Proposition[ 
ensures that a n — > 00 as n — > 00 and a similar consideration gives b n — > 0. Consequently the 
angle between vectors v n and tends to zero in the limit, with the former strictly in C and the 
latter by assumption not. The implication is that both approach the boundary of C, given by 
(dx,dy) £ span{(l,0)}Uspan{(0, 1)}. 

However cannot approach (1,0), for consider some iteration z = H m (z) so that H (z) = 
F{z). Then 



DH = DF = 



1 

D x f+{z) D 2 f+(z) 



By continuity of the linear map DF any neighbourhood U of (1, 0) is mapped into a neighbourhood 
V of (1, £>i/+(z)). Because Dif + (z) is bounded uniformly away from zero (as determined in the 
previous section), one can find U so that V C C. So e" cannot approach (1,0) as in doing so it 
is inevitably mapped into C, contradicting the assumption. Analogously one can show that 
cannot approach (0, 1) by considering some iteration z = H m (z) so that H (z) — G (z). This gives 
the required contradiction. □ 

By similar arguments E s (z) £ C(z) and the length of H~ m (7 s (z')) diverges to infinity as 
m — ► 00 . The orientations of (un)stable subspaces have the immediate consequence that gradients 
of (un)stable manifolds are similarly aligned on the manifold itself. The remainder of our proof is 
essentially geometric and is simplified by the introduction of a covering space for R. We do this 
in two stages. 

Let —R C T 2 denote those points (x, y) C T 2 so that (~x, — y) £ R. Notice that R n -R = 
and let Ri — R U —R. Define p' : Ri — > R by (x, y) 1— * (x, y) if (x, y) £ R and (x, y) 1— > (—x, —y) 



A Bernoulli linked-twist map in the plane 



12 



otherwise. Then (f? 1; p') is a covering space (a double cover, in fact) of R. The derivatives of p' 
and its possible inverses each preserve the cones C and C. 

Let 7r : M 2 — > T 2 be the natural projection which takes each coordinate modulo 2tt and let 

R 2 = {(u,v) G R 2 : ir(u,v) G i?i} C M 2 . 

i?2 has the form of a lattice, constructed by fitting together an infinite number of copies of R\ and 
is illustrated in Figured] Let p" : R2 — > R\ be the projection which takes each coordinate modulo 
S 1 . Then (i?2,p") gives a covering space for R\, The derivatives of p" and its (local) inverses 
preserve C and C. So (i?2,P = p' °p") is a covering space for i?, and Dp, Dp -1 preserve C and C. 
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Figure 8. A portion of the manifold R2 C K 2 . Together with the map p : R2 —> R this gives 
a covering space for _R. In red is a typical piece of some image of a local unstable manifold for 
some 2 S R. The gradient at all times is in C. Analogously in blue is a typical piece of some 
pre-image of a local stable manifold for z' G R. Its gradient is in C. If each is sufficiently long 
then they must intersect. 

It is now elementary to show that (|7.2[) is satisfied, for if we consider any sufficiently long 
H m (7"(z)) and any sufficiently long H ~ n (7^(21')) we can always lift them to i?2 in such a way that 
they intersect. Figure [5] illustrates an example. The image with respect to p of the intersection 
point is an intersection point in R. This completes the proof of Theorem ll.il 

8. Concluding remarks 

Our method necessitates a strong restriction on the sizes of the annuli on which the linked-twist 
map is defined. The restriction is used in proving the DiJ-invariance of the tangent cone C 
(Proposition 16 . lj ) where, for example, it was required to show that 
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>0, 



r 27r 

D x ijj (x,^~ 1 (x,y) + c(x - r )) + D 2 ip (x,tf}~ 1 (x,y) + c(x - r )) D\ij)~ {x, y) H 

L ^1 — ^0 

for each pair (x,y) <G [?"o,7t] x [0, 7r] . Analytically determining tight estimates on the left-hand 
side is very difficult (even plotting it using computer algebra software requires a non-trivial effort 
because of the different forms taken by tp and its inverse). By comparison our approach of bounding 
each of D\i()(-, •), D2ip(-, •) and ■) individually is rather crude. Although the lower bound 

of for D±ip is optimal, none of the other bounds established are. ft is remarkable that the 
DiJ-invariance of C may be established for any choice of annuli using this approach and perhaps 
indicative that suitable bounds in fact hold for a much wider choice of annulus size. One obvious 
way to resolve this is to partition the domain of the functions so that tighter bounds can be 
established element-wise; the aforementioned computer plots might suggest a sensible partition. 

We consider how far one might proceed in this manner. It is natural to wonder whether 
one can determine a set of values for r$ and r% from which the Bernoulli property follows, and a 
complementary set on which it is shown not to occur. Numerical simulations of Sturman et al. 
(2006) would suggest that the present method is insufficient for this task for the following reason. 
Essential to our ability to construct the new coordinates is that £± are disjoint, so that the shears 
of the respective twist maps act transversally for each point in £. The simulations suggest that 
such transversality is not a necessary condition for good mixing. It remains an interesting open 



question as to whether transversality and Wojtkowski's condition (|3 - 2[) are sufficient for Bernoulli 



Acknowledgments 

Whilst conducting this work JS was supported by EPSRC and SRW by ONR Grant No. N00014-01- 
1-0769 and EPSRC Grant EP/C515862/1. The authors are grateful to Rob Sturman for a careful 
reading of the text, and to Holger Waalkens and Jens Marklof for many helpful suggestions. 



References 



Bowen R 1978 number 35 in 'Proc. CBMS Regional Conf. Math. Ser.' Amer. Math. Soc. Providence. 

Braun M 1981 SI AM J. Math. Anal. 12(4), 630-638. 

Chernov N & Haskell C 1996 Erg. Th. Dyn. Syst. 16(1), 19-44. 

Devaney R L 1978 Proc. Amer. Math. Soc. 170, 71(2). 

Devaney R & Nitecki Z 1979 Comm. Math. Phys. 67(2), 137-146. 

Katok A, Strelcyn J M, Ledrappier F & Przytycki F 1986 Invariant Manifolds, Entropy and Billards; Smooth Maps 

with Singularities Vol. 1222 of Lecture Notes in Mathematics Springer- Verlag, Berlin, New York. 
Liverani C & Wojtkowski M 1995 Dynamics Reported (New Series) 4, 130-202. 

Moser J 1973 Stable and Random Motions in Dynamical Systems Princeton University Press Princeton. 

Ottino J M 1989 The Kinematics of Mixing: Stretching, Chaos, and Transport Cambridge University Press 

Cambridge, England. Reprinted 2004. 
Ottino J M & Wiggins S 2004 Science 305, 485-486. 
Pesin Y B 1977 Russ. Math. Surveys 32, 55-114. 

Przytycki F 1981 Linked twist mappings: Ergodicity. Preprint, IHES. 
Przytycki F 1986 Studia Math. 83, 1-18. 

Springham J 2008 Ergodic properties of linked-twist maps PhD thesis University of Bristol. ArXiv:0812.0899vl. 
Sturman R, Ottino J M & Wiggins S 2006 The mathematical foundations of mixing Cambridge University Press 
Cambridge. 

Thurston W P 1988 Bull. Amer. Math. Soc. 19(2), 417-431. 

Wiggins S & Ottino J M 2004 Phil. Trans. Roy. Soc 362(1818), 937-970. 

Wojtkowski M 1980 in 'Nonlinear dynamics (Internat. Conf., New York, 1979)' Vol. 357 of Ann. New York Acad. 
Sci. pp. 65-76. 



